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Abstract 

Recent observation shows that the Higgs mass is at around 125 GeV while the prediction of the 

minimal supersymmetric standard model is below 120 GeV for stop mass lighter than 2 TeV unless 

the top squark has a maximal mixing. We consider the right-handed neutrino supermultiplets as 

messengers in addition to the usual gauge mediation to obtain sizeable tridinear soft parameters 

At needed for the maximal stop mixing. Neutrino messengers can explain the observed Higgs mass 

r*j ■ for stop mass around 1 TeV. Neutrino assistance can also generate charged lepton flavor violation 

P-I 

including [i — > ej as a possible signature of the neutrino messengers. We consider S4 discrete flavor 



model and show the relation of the charged lepton flavor violation, #13 of neutrino oscillation and 
muon g — 2. 
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I. INTRODUCTION 



The observation of the Standard Model Higgs-like new boson with mass at around 125 
GeV changes the current understanding of new physics at the weak scale. The minimal 
supersymmetric standard model (MSSM) can explain 125 GeV with a relatively light stop of 
1 to 2 TeV in the context of maximal stop mixing. From the model building point of view, it 
is quite difficult to realise the maximal stop mixing scenario starting from ultraviolet (UV) 
theory. In minimal gauge mediation (MGM) (3 6 1 , soft tri-linear A term is not generated 
at the messenger scale and the radiatively generated A term at the weak scale is not large 
enough to realise the maximal stop mixing. As a result, colored superpartners should be as 
heavy as 5 to 10 TeV to explain 125 GeV mass of the Higgs boson 7, sj]. Therefore, the 
explanation of 125 GeV Higgs boson mass needs an extra help in minimal gauge mediation. 
Next to the minimal supersymmetric standard model (NMSSM) can use the extra contri- 
bution from the Yukawa F-term of the singlet. For instance, look at Extra vector-like 



fermions are added in minimal gauge mediation 10Ml2|. Direct coupling of visible sector 
fields with messengers can help. Higgs-messenger mixing jl^, 1^] and matter- messenger 



mixing 
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19j can generate Yukawa mediated contribution including A term at the mes- 



senger scale. However, at the same time the virtue of gauge mediation is gone and they 
would spoil nice flavor preserving spectrum and can possibly cause the flavor problem at 
the weak scale. For the Higgs-messenger mixing, A/m? problem (13I. \i\ which is analogous 
to \1jB\1 problem can arise and the electroweak symmetry breaking is difficult to achieve if 
the mixing coupling is large. General gauge mediation 20| can avoid this problem by using 
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22] 



the mechanism of radiatively generated maximal stop mixing 

In this paper we consider the right-handed neutrino supermultiplets as the messengers of 
supersymmetry breaking in addition to the messengers charged under the Standard Model 



(SM) gauge group, e.g. 5 and 5 of £77(5). The setup is motivated from [23| which provides 



a solution to \i problem in gauge mediation (more precisely \ijB\i problem) |24l . |25|. For 
the solution in [23[ to work, the messenger scale should be higher than the Peccei-Quinn 
breaking scale, 10 9 ~ 10 11 GeV. For the gauge mediation to be the dominant contribution 
compared to the Planck suppressed higher dimensional contribution, the messenger scale 
should be lower than 10 15 GeV. Therefore, the See-Saw scale with order one neutrino Dirac 



Yukawa couplings, 10 13 ~ 10 14 GeV, is well motivated as the messenger scale if we accept 



23J as a solution of the /x problem in gauge mediation. If the messenger scale is at around 
the See-Saw scale, the natural question is why the right-handed neutrino supermultiplets 
do not serve as messengers of supersymmetry breaking. Apparently there is no harm to 
couple the right-handed neutrino superfields directly to the messengers. Majorana mass of 
the right-handed neutrino and the messenger mass of ordinary 5 and 5 might have the same 
origin in this case. In summary, the minimal set of messengers are 5, 5 and 1. This is 
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different from previous studies relating gauge mediation and See-Saw mechanism 
They employ particles relevant to See-Saw mechanism as messengers, and these particles 
are also charged under the SM gauge group, which can be seen in the Type-II or Type-Ill 
See-Saw. Therefore, gauge mediation and neutrino Dirac Yukawa mediation have a common 
messenger. In our case, in contrast, neutrino Dirac Yukawa messenger is the right-handed 
neutrino, the SM singlet. 

If the right-handed neutrinos couple to the supersymmetry breaking field, neutrino Dirac 
Yukawa coupling generates the A term and soft scalar mass of lepton doublet and up-type 
Higgs at the See-Saw scale after integrating out the neutrino messengers. 125 GeV Higgs 
mass can be explained with stop lighter than 2 TeV in this setup. At the same time the 
stop mass gets an extra Yukawa mediation and maximal stop mixing can be easily realised. 

As there is a neutrino Dirac Yukawa contribution to the soft parameters in addition to the 



ordinary gauge mediation, interesting new physics signature is expected. T 
;he charged lepton flavor violation is different from that in mSUGRA 



re mechanism of 



291 or SUSY GUT 



30j in which the origin is the running of soft parameters above the See-Saw scale. Though 
the origin is different, the spectrum looks similar. The crucial difference is that here the 
flavor violation appearing in lepton doublet soft scalar mass is IQtt 2 bigger than the one 
in mSUGRA. Therefore, the naive expectation is that order one neutrino Dirac Yukawa 
coupling would be incompatible with the current bounds of various charged lepton flavor 
violation constraints including /i — > ej. 

The computation of the charged lepton flavor violation needs a complete flavor model. 
Current observation of the charged lepton mass and lepton mixing matrix (PMNS) can be 
explained in a consistent way with the neutrino Dirac Yukawa matrix which is proportional 



to the identity matrix. This is not an ad hoc assumption but can be exp 
context of non-Abelian discrete flavor symmetry, e.g., tribimaximal PMNS 



ained in the 



311 ] from S4. 



Therefore, order one neutrino Dirac Yukawa coupling can generate order one A term at 



the messenger scale and at the same time can be consistent with the charged lepton flavor 
violation constraints as long as it is proportional to the identity matrix. 

£4 flavor symmetry is the most natural and/or simple if #13 = as the tribimaximal 
mixing can be nicely realised. However, small but sizeable #13 (sin 6^3 ~ 0.15) can be 
accommodated with the extra complication 32M35| . If the origin of #13 is the modification of 
Majorana mass of the right-handed neutrino, there would be no off-diagonal element in the 
lepton doublet soft scalar masses as the neutrino Dirac Yukawa would be still proportional 
to the identity. In this case the model is free from the cLFV constraints. Nevertheless, the 
sparticle spectrum needed to explain the observed Higgs mass is heavy enough such that it 
is hard to explain the muon anomalous magnetic moment at the same time. If #13 is due to 
the deviation of the neutrino Dirac Yukawa matrix from the identity, sizeable charged lepton 
flavor violation is expected. We compute the charged lepton flavor violating processes in 
both cases and show that interesting parameter space exists if #13 is a combination of two 
contributions from neutrino Dirac Yukawa and Majorana mass matrix. 

The contents of the paper is following. In section 2, we explain the setup for neutrino 
assisted gauge mediation in which the right-handed neutrino is added as messengers in 
addition to the ordinary SM charged messengers. Also we discuss the implication for the 
Higgs mass. In section 3, we explain our £4 flavor model as a representative example to 
discuss possible phenomenological implication. In section 4, we discuss charged lepton flavor 
violation in connection with muon anomalous magnetic moment, the neutrino mixing angle 
#13 and the Higgs mass. Then we conclude. 



II. NEUTRINO ASSISTED GAUGE MEDIATION AND THE HIGGS MASS 



A. Soft terms generated from right-handed neutrino messengers 



The extremely small masses of neutrinos can be explained through the See-Saw mechanism 36 



401 ] . in which lepton number is violated at around the Grand Unified Theory (GUT) scale. 
In this paper, we consider the simplest model, type-I See-Saw. For this, we extend the 
MSSM superpotential by including right-handed Majorana neutrinos, 

a. zjb (\r \ n n« Tjb l\r \ fp to, rjb , (\r \ at r a jjb 



W = e ab 



{Yu^UiQTjK ~ (Yn^D^Hl - (Y E )ijE i L a j H d + (Y^jN^H* 

1 - 
2" 



+ -MSNiNj, 



where e a b is a totally antisymmetric tensor with 612 = 1. The superfields in the superpotential 
represent right-handed neutrino-sneutrino pairs, in addition to the SM particles and their 
superpartners. They have the following SM gauge group SU(3) c xSU(2) £ xU(l)y quantum 
numbers: 

Q:(3,2,~), t7 : (3, 1,-|), D : (3, 1 i) 
o 3 3 

L:(l,2,-i), £: (1,1,1), iV : (1,0,0) (2) 

ff„:(l,2,i), ff d :(l,2,~). 

Relative minus signs of Yukawa terms are given to make the sign of terms responsible for 
the fermion Dirac masses to be the same. 

The relevant soft supersymmetry (SUSY) breaking terms are given by 

£soft = - Hv)}NjN, - {m\)\m 3 - m\EiR u 

^(B N M) ij NiNj + {A v ) l3 U l QiH u - (A D ) ij D i Q j H d - {A E ) lJ E l VH A + BfiH u H d + h.c. 

(3) 

We consider two origins of soft terms. The first one is gauge mediation. In the gauge 
mediation, sfermions obtain soft masses given by jsj 

F \2 



m 2 



h^(ii)«.i:y (4) 

a i 

at the messenger scales M i; where C a is the quadratic Casmir ^ a T Q T a of the sfermion 
representation IZi under the corresponding gauge group labeled by a, which is given by 
(N 2 - l)/(2iV) for SU(N) and Y 2 for U(l) y , T a is defined by TrT ^ = T a (^)5 a ^, and 
f(xi) is the loop function of Xi = F/M 2 which is close to one for small Xi. On the other 
hand, the 2-loop tri-linear A term is very small and can be neglected at the messenger scale. 

For the second origin of soft terms, we introduce a SUSY breaking spurion X which 
couples to the right-handed neutrinos. Majorana mass of the right-handed neutrino comes 
from the scalar vacuum expectation value (VEV) of the SUSY breaking spurion X, 

W D XXNN. (5) 

Then N acts as the messengers of supersymmetry breaking, and the neutrino Dirac Yukawa 
coupling, 

W D Y U NLH U , (6) 
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is interpreted as the direct mixing term among the messengers, Higgs and matter (leptons). 

The SUSY breaking effects at the See-Saw scale Mm = A(X) is studied in |4l|. When 
right-handed neutrinos couple to the SUSY breaking sector, Majorana mass matrix is ana- 



lytically continued to be — > (1 + 9 2 B^)M^, as in the case of gauge mediation 42N44 1 . 
Here, we assume that the flavor structure of the right-handed neutrinos is fully determined 
by M^, so Bn = Fx/X is a constant. 

Then, SUSY breaking is transferred to the visible sector through the neutrino Dirac 
Yukawa interaction. Wave function renormalization from the interaction with right-handed 
neutrinos is given by 

where 



16vr 2 V A 2 



A* = [Z~ 1/2 ] t *nZ2 1/2 ZhI /2 , M r = [Z-^fM^ 2 , (8) 



then analytically continued Majorana masses give the soft masses. From field redefinitions 

(9) 



L^(l- 5 -^)(l- 6 2 6Z L \ B2 )L 



H u ^(l- 6 -^)(l-6 2 5Z H MH ul 
supersymmetric kinetic terms can be written in the simple form, 

$ f (i + <*z*)$ + eWsz^i^)® (io) 

then we can read off the one-loop corrections to the soft masses 

5m\ = -5Z L \ e2 g 2 and 5m 2 Hu = -SZ Hu \ 92 §2. (11) 

In the expression, is just a constant, not a matrix. So m(Mjy-Mjv) in the wave func- 
tion renormalization is separated into holomorphic and anti-holomorphic parts, respectively. 
Since 9 2 9 2 term is not generated, we do not have one-loop soft masses. 

Hence, as in minimal gauge mediation, soft masses are generated at two loop level. In 



13|, it was shown that soft scalar masses of the fields which directly couple to messengers 
and those which do not are different. In our model, the slepton L and the up- type Higgs 
H u couple to messengers iV directly to give soft terms, 
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5m 2 T 



B 



N 



(An) 



3g 2 2 - -gl)YX + 3^7,^7, 



5m 2 



D 



N 



ATr[Y v YWY v 



3gl+-gl)Tr[Y„Yj] 



;i2) 



(4tt) 4 L 

On the other hand, Q and U obtain two-loop soft scalar masses through the wave function 
renormalization of H u and the corrections are given by 



5mi 



dirty 



D 



N 



(4tt)< 

B N 



(4tt)< 



(l\[Y v Y^Y v 
-MY^YuYt 



(13) 



while the soft masses of E and Hd come out of the wave function renormalization of L and 
the corrections are given by 



5m 2 



E 



Sm H d 



b2n y,y ] y y ] 

(4tt) 4 v e 



(14) 



By replacing Y E -» Y E (1 + 8A E ), Y v ->■ Y L1 {1 + 5A V ), and Y D ->> Yd(1 + <L4 D ), we have 
following soft terms at one loop level, 

8A E — —6Z L \ 9 2, 5Au = -I5Z Hu \e 2 , 
5A D = 0, 5B = -5Z Hu \ 6 2. 

Unlike gauge mediation, right-handed neutrino mediation generates one-loop A— terms, 

B N 



(15) 



A, 



.yty 

, 1 V 1 v 



167T 2 

A v = -TrA E x I 3x3 (16) 
B = TrA E . 

While gauge mediation contributions are flavor universal, See-Saw Yukawa mediation 
is flavor dependent and one of the virtue of the gauge mediaion would disappear. In the 
absence of See-Saw Yukawa mediation, cLFV can appear when the messenger scale is higher 
than the right-handed neutrino Majorana mass scale. See-Saw Yukawa contributes to the 
slepton soft mass through the renormalization group equation (RGE) , 

1 



d q d o 
lx—m T = n—rrir 
dfx a fx 



MGM 16"7T 2 



m 



\YX + Y}Y v m\) + 2{Y}m 2 N Y v 



m 



Y}Y U + At A 



(17) 
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which should be restricted by cLFV constraints 45[. Here A v = A V Y V is used. Since m\ is 
two-loop generated, cLFV effects are further loop suppressed (at three loop level). Unlike 
mSUGRA, this effect is known to be small in gauge mediation as the messenger scale is at 
most comparable to the See-Saw scale and the running can be made in a very short interval. 
This is not the cLFV that we are interested in. 

In neutrino assisted gauge mediation, neutrino Dirac Yukawa couplings can introduce 
two-loop generated cLFV effects on m\ as a result of gauge- Yukawa or Yukawa mediation, 



5m 



2 N 



(47T) 



(Tr^yj] + 3Tr[Y^] - 3g 2 2 - \gf)Y^ v Y v + SY^X] (18) 



in the charged lepton mass basis. If the two loop generated slepton mass squared has a 
nonzero off-diagonal element, it would generate cLFV. Parametrically, this effect is much 
larger than the expected cLFV in mSUGRA or similar scenarios in which the effect comes 
from the running above the See-Saw scale. We simply assume that both messengers 5, 5 and 
1 have the same masses at the See-Saw scale. In principle these two masses can be different 
and cLFV can arise if the singlet messenger is lighter than 5,5. However, this effect is loop 
suppressed compared to the Yukawa mediation we would not consider it in this paper. 

Further discussion on cLFV is possible only when there is an explicit flavor model pro- 
viding the neutrino Dirac Yukawa and charged lepton Yukawa matrices. As a simple and 
illustrative example of the explicit model, we consider S4 flavor symmetry in Sec. Ill II It will 
be shown that various types of See-Saw Yukawa Y v would predict different sizes of effects 
on cLFV. Before moving onto the flavor discussion, let us consider the implication on the 
Higgs mass first. 



B. Higgs mass and superparticle spectrum 

Minimal gauge mediation does not generate A t at one loop and the weak scale A t is 
radiatively generated by the gluino loop. However, the same gluino contribution appears in 
stop soft scalar mass and the relative ratio of | A t \ and mi can not be large. On the other hand, 
the physical light CP even Higgs mass in the MSSM is affected by X t = (A t — \ij tan /3)/m 4 - 
and X t ~ 2 (or \/6 more precisely) gives the maximum finite threshold correction as shown 
in Fig. [TJ 

One way to make \X t \ > 1 at the weak scale is to start from tachyonic stop boundary 
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m,,(GeV) 




FIG. 1. Higgs mass with respect to Xt for tan/3 = 10, stop mass Mj ~ 2TeV. 



condition 2l| which is explicitly realised in gauge messenger model 22|. However, this option 
is not available in minimal gauge mediation. The other possibility is to couple messengers 
directly to the visible sector fields such that large negative A term can be generated at the 
messenger scale. If A term is positive, the gluino contribution from the running cancels the 
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also has been considered 



A term at the messenger scale. Matter-messenger mixing 
recently. Messenger-matter-matter Yukawa coupling would generate the needed A t term 
at the messenger scale. However, the full Yukawa couplings are written as 3 x 3 matrices 
and why all other dangerous Yukawa couplings between matters and messengers are absent 
except 33 component remains to be a puzzle. One way out is to consider Higgs-messenger 
mass mixing 44[ and to generate all the wanted Yukawa couplings between matter and 
messengers from ordinary Yukawa couplings of matter with Higgs. There would be a direct 
modification of squark spectrum if squark couples directly to the messenger. 

Higgs-messenger mixing through Higgs- messenger- messenger coupling or Higgs-Higgs- 
messenger coupling has been considered in 13J, LL4{. In this case, we often encounter 
A/m 2 problem. To understand this, it is worth to emphasize that, the two-loop soft mass 
squared of the Higgs field H u which has direct coupling to messenger $ has a structure of 
m 2 H ~ cA 4 — c'X 2 g 2 where A is a coupling constant of Higgs and messenger fields and g is the 
gauge coupling(s). On the other hand, the two-loop soft mass squared of fields Q, U which 
does not have a direct coupling with messenger has a form of ^ ~ —c^\ 2 y 2 . This fact 
was extensively studied in [13j. For sufficiently large A, large one-loop A terms are generated. 
At the same time, m 2 H becomes positive so the soft mass of H u can be much larger than 
that in the pure gauge mediation. Moreover, the soft mass of Q3, U$ can be much smaller. 
If the Higgs H u superfield directly couples to messengers whereas the top superfields do not, 



M, 



,rannmg 



164 



162 



Stdp Tacrjydnic; • 



160- 



158 



0.66 



M, 



f,runrung 



164 



162 



160 



15b 



0.68 



0.70 



■ .*■•*■ ■" 



No ewsb : -v. 

0.72 ' " 0.74 i_ 0.763 ; v 



■ ':Stop'Tachyohic : 



..No.BWSB; ; 



FIG. 2. Phase diagrams indicating stop tachyonic and no EWSB region for tan/3 = 10, tan /3 = 30, 
respectively. is set to be 5 x 10 5 GeV. 

relatively light stop in natural SUSY can be easily obtained as we can have the small stop 
soft mass from the effect explained above and the large LR mixing from large A term. A/m? 
problem appears in H u soft terms such that large A term at the same time generate large 
m 2 Hu at the messenger scale and it can make the electroweak symmetry breaking difficult. 
It is analogous to the famous n/Bfi problem in gauge mediation. To avoid this but to allow 
the large A for maximal mixing, large —C2\ 2 g 2 part in m 2 Hu is required. This can be achieved 
by introducing new gauge bosons or making strong interaction involved On the other 
hand, one loop, negative contribution to rn 2 Hu can be considered if the messenger scale is low 
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FIG. 3. Higgs mass as a function of y v for = 5 x 10 5 GeV, p = 0.1. Higgs mass can be achieved 
with the help of Yukawa mediation for large tan f3 region. At y v ~ 0.7, stop mass is approximately 
ITeV. 

as analysed in detail in jl^ ]. 

Neutrino assisted gauge mediation uses the Yukawa coupling among messengers (neutri- 
nos), Higgs and lepton doublets. As a result, Higgs and lepton doublet soft scalar masses 
get extra contribution from Yukawa mediation. The same A/m 2 problem applies here and 
neutrino Dirac Yukawa coupling can not be taken to be a large value for successful elec- 
troweak symmetry breaking in principle. On the other hand, too large m 2 H and too large 
A term may drive stop tachyonic through renormalization group running with top Yukawa. 
The problem becomes worse as the stop soft scalar mass squared at the messenger scale 
gets a negative contribution from Yukawa mediation. The situation is shown in Fig. [2j For 
the running mass of the top quark 160 GeV (the central value), the tachyonic stop appears 
before the real A/m 2 problem prevents the successful electroweak symmetry breaking as we 
increase y v . The crucial difference compared to the previous work in which A/m 2 problem 
is emphasized [3, [ljj] comes from the number of messengers. In neutrino assisted gauge 
mediation, the number of messengers is three, N = 3. The y 2 contribution is effectively 
replaced by Ny 2 . Large N effectively reduces the A/m 2 problem by 1/N . At the same time 
smaller y u can provide the same impact with the aid ot N > \. If tachyonic stop appears as 
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y v gets larger, it would be easy to realise the maximal stop mixing by making the stop soft 
scalar mass sufficiently small. 

3500 ' 
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FIG. 4. M~ t as a function of y v for B N = 5 x 10 5 GeV, p = 0.1. tan (3 = 10. 
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FIG. 5. X t /M- t as a function of y v for B N = 5 x 10 5 GeV, p = 0.1, tan /3 = 10. 

Fig. [3] shows the contribution assisted by neutrino messengers, compared to the minimal 
gauge mediation which corresponds to y v = with stop mass at around 1 TeV. In the 
minimal gauge mediation, the Higgs mass is computed to be at around 121 ~ 122 GeV for 
tan — 10 ~ 30. For y v = 0.7, the Higgs mass can be as large as 125 ~ 126 GeV. 4 to 5 
GeV gain in the Higgs mass is obtained in neutrino assisted gauge mediation. The gain does 
not look impressive but has an impact on allowed superparticle spectrum. In the absence of 
At at the messenger scale as is the case in minimal gauge mediation, this extra 5 to 6 GeV 
can be achieved by making the logarithmic contribution large and the stop mass should be 
as heavy as 5 to 10 TeV rather than 2 TeV. 
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FIG. 6. At = A t Y t as a function of y v for tan/3 = 10, = 5 x 10 5 GeV. Without Yukawa 
mediation, one can obtain At ~ — 2700 GeV at the weak scale by RG running effects. With help 
of neutrino mediation at the messenger scale, one can obtain At ~ —4500 GeV at weak scale. This 
drives more stop mixing, which helps 125 GeV Higgs mass. 

Note that the plot stops at y v = 0.72. Neutrino assisted gauge mediation is classified 
as Higgs-messenger mixing scenario as the right-handed neutrino is the messenger and the 
neutrino Dirac Yukawa coupling connects Higgs, lepton doublet and the messenger (right- 
handed neutrino). The stop soft scalar mass squared gets smaller and becomes tachyonic as 
the neutrino Dirac Yukawa coupling is increased as in Fig. HI The logarithmic correction 
to the Higgs mass also rapidly drops beyond y v ~ 0.7 as the stop mass becomes too light 
(and becomes tachyonic) as is shown in Fig. [3j The maximal mixing is realised around this 



point, as shown in Fig. |5j T 



to the critical point as in 46]. 



lis also corresponds to the corner of the parameter space next 



Fig. [6] compares A t in the minimal gauge mediation and the neutrino assisted gauge 
mediation both at the messenger scale and the weak scale. Note that A t by itself is enhanced 
by 1.5 at the weak scale with the help of messenger scale A t . 

Fig. [7] shows the relation between Bn and the Higgs mass. The neutrino Dirac Yukawa 
coupling y v is chosen to be close to 0.72 which can maximize the Higgs mass for given B^. 

In summary, the minimal gauge mediation needs stop mass at around 5 to 10 TeV to 
raise the Higgs mass up to 125 GeV. If the right-handed neutrinos are the messengers of the 
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FIG. 7. Maximum values of Higgs mass as a function of Bjy. For tan (3 = 10, at least = 360 TeV 
is required to obtain 125 GeV Higgs mass. For tan/3 = 30, = 300 TeV is required. At two 
points (tan/3 = 10, B^ = 360 TeV), (tan/3 = 30, Sat = 300 TeV), we display sparticle spectrums 
in Table ITTT1 Also spectrums with 123 GeV Higgs mass are given for (tan/3 = 10,2? at = 240 TeV), 
(tan/3 = 30, Bn = 200 TeV). Sparticle spectrums are displayed in Table ITVl 

super symmetry breaking, so called 'neutrino assisted gauge mediation', we can explain 125 
GeV Higgs mass with lighter than 2 TeV stop mass. 

III. FLAVOR MODEL 

In this section, we consider models which can explain neutrino oscillations successfully. 
Since the SUSY breaking mediation through the neutrino Dirac Yukawa coupling is flavor 
dependent in general, sizable cLFV could be generated. To avoid this, the neutrino Dirac 
Yukawa coupling is set to be proportional to the idenitity. In the right-handed neutrino 
mass basis, it would be proportional to the unitary matrix so soft mass m 2 L , which de- 
pends on the combination YjY v is flavor universal. It is easily achieved by employing the 
non-abelian discrete symmetry for the tri-bi maximal mixing of the PMNS matrix. Since 
the tri-bi maximal mixing should be modified to make #13 nonzero, as reported by several 
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observations 47H51|. small corrections should be added. When the neutrino Dirac Yukawa 
coupling has such corrections, such that it has a deviation from identity, cLFV is generated. 
We look for several ways to suppress cLFV, at least under the experimental bound. 
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TABLE I. Charge assignments under S4 x Z4 x U(1) L x U(l)^ for leptons, flavons, Higgs, and 
SUSY breaking spurions. 



To make the PMNS matrix tri-bi maximal, we use S± discrete symmetry, since it is closely 
related to the permutation structure of Yukawa couplings. Other discrete symmetries, such 
as A4, the even permutation of the S4 could be used. The main difference is that the first and 
the second generation of the right-handed leptons belong to 2 dimensional representation 
in 5*4 while they correspond to different one dimensional representations, 1', 1" in A4. In 
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55|, the structure we use is obtained from A4 symmetry and discussion on the deviation 
from the tri-bi maximal mixing is in parallel. The S4 symmetry model building is reviewed 
in js6|. In Appendix A, we summarised representations and tensor products of S4 group. 

For quark sector, the CKM matrix is close to the identity. Deviation from the identity 
has a hierarchy structure parametrized by some powers of the Cabibbo angle, A = sinOc- 
On the other hand, the PMNS matrix, mixing matrix in the lepton sector has large mixing 
angles. Even the smallest mixing angle, #13 is in the order of A. To explain this, it is natural 
to assume that u— and d— quark sectors have almost the same structure under the discrete 
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flavor symmetry whereas the charged lepton and the right-handed neutrino sectors do not. 
This picture can be realised by introducing appropriate 'Savons' charged under discrete 
symmetry group and more symmetries can be introduced to forbid useless couplings. Here, 
we consider the symmetry group S4 x Z 4 x U(1) L , where U(1)l represents a lepton number, 
which may be discretized. In this paper, we consider superpotential for See-Saw mechanism 
with flavons $ and Xi 

W = + hijNiLjHu + U 3lj XN iX N v (19) 

where i, j = 1, 2, 3 are the generation indices and X is a SUSY breaking spurion. For this, 
S4, Z4, U(1)l and U(l) R-symmetry quantum numbers are given in Table [H 

The charged lepton Yukawa couplings can be constructed from EQLHd, the neutrino 
Dirac Yukawa coupling can be constructed from NL, and the Majorana mass of the heavy 
neutrinos can be constructed from XN\N . On the other hand, <J> 2 , y 2 , and $x cannot 
couple to the combinations ELH^, NL, and XNN to make singlets. Note that U(l)^ is 
introduced to forbid unwanted coupling NyN, which makes -Bat in a matrix form, not a 
constant. 

The discrete symmetry quantum number can be extended to the quark sector, such as 
Q : (3, 1, 1, 1, 1), U : (2 + 1, 2, 0, 0), and D : (2 + 1, 2, 0, 0) under S 4 x Z 4 x U(l) L xU(l) i? . 
The flavons $ : (3 + 3', 1, 0) make the singlet combinations U§QH U + D^QH^ and Yukawa 
couplings Yjj and Yd have the same form as the charged lepton Yukawa coupling. They are 
diagonalized by the same unitary matrix so CKM matrix is the identity in the leading order. 
If another type of flavon couples to either of up and down quark sectors to give subleading 
corrections of order A, it would explain the Cabibbo angle. 

Lepton Li is in the 3 and Ej is in the 1+2 representations, in which (£1)1 + (E 2 ,E 3 )2- 
Also there are the SM singlet flavons $3, and $3' in the 3, and 3' representations. We do 
not provide a complete vacuum alignment in this setup. Instead in Appendix B, we show 
a few simple examples in which the aligned vacuum is realised. If, for instance, VEVs are 
arranged to be ($3) = w 2 (l,l,l), and ($3') = ^3(1, 1,1), we have the following Yukawa 
structure 



( c c c ^ 



a au au 2 



b bu 2 bu 



(20) 



1(3 



where a = (At>2 + X'vs)/A, 6 = (XV2 — X'v^j/A, c = A'^/A, and A, A', A" are coupling 
constants of E2L3&3, E 2 Ls^ S ', and ExL^s, respectively. In this case, Y e Ye has the form 
of 



/ 



Y E Y E 



\X, 



c 2 + a 2 u + b 2 uj 2 c 2 + b 2 oo + a 2 ui 2 



\ 



\ 



c 2 + a 2 u 2 + b 2 u c 

c 2 + b 2 co 2 + a 2 u c 2 + a 2 J 2 + b 2 cu a 2 + b 2 + c 2 
which will be diagonalized to |A^| 2 ((e 3 ) 2 , (e) 2 , 1) by the unitary matrix, 

(11 A 



b 2 + c 2 c 2 + a 2 u + b 2 oj 2 



(21) 



1 w 2 w 
1 cu u 2 



(22) 



Here we use e ~ m At /m T as the order parameter. Then, c = e 3 , a = e and 6=1. 

On the other hand, let heavy neutrinos iVj be in the triplet 3. $3 and $3/ cannot couple 
to the combination LjiVj by Z 4 and U(l)^ symmetries as well as the SM gauge symmetry. 
Since the combination L\N\ + L2N2 + L3N3 is a singlet, we naturally have the neutrino 
Dirac Yukawa coupling Y v proportional to the identity. Finally, XNiNj has again the form 
of 3 + 3' + 1 + 2. $'s cannot couple to it while singlet xi an d triplet xs i n the singlet and 
triplet representation can do, so we have the following Majorana mass term: 



^ Wi W2 ^ 



M 



N 



wi 
w 2 Wi 



(23) 



where (xi) = n>i and (xs) = w 2 (0, 1,0), respectively. Therefore, the neutrino mass matrix 
My = — ^YjfM^Yy is diagonalized by 



/JL -JL\ 
1 

V75 71 / 



(24) 



so we obtain the PMNS matrix in the tri-bi maximal mixing 

/ 



1 

1 



„— i5w/6 1 



V 



V2 

,,2 1 ,,2 1 p i57r/6 1 

- u 75 w 7s e 73 / 



(25) 
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In this construction, S4 triplet flavons have VEVs in the direction of (1,1,1) or (0,1,0). 
These directions are easily stabilized compared to other directions, such as (1, 1, 0), as argued 
in Appendix B. Note that Y v proportional to the identity does not give rise to LFV. In Eq. 
f TTZj) . we see m 2 L from the neutrino Dirac Yukawa mediation is flavor universal. In the right- 
handed neutrino and the charged lepton mass basis, Y v moves to VjY v and m 2 L moves to 
{Viym 2 L V[. As a result, PMNS matrix is mulitplied and will change m 2 L matrix. However, 
if the neutrino Dirac Yukawa matrix is proportional to the identity matrix, the property of 
VW ] = WV ] = I cancels out such effects. 

There are various ways to put corrections to make non-zero 8i 3 . Moreover, corrected 
neutrino mass matrix should be consistent with the measurements of 612, 623 as well as 



neutrino mass squared differences, Am 



2 

sol 



TTln 



m\ and |Am 



2 

atm I 



\nio 



m 2 .]. Since the 



overall neutrino mass scale is not known, the important quantity is the ratio of neutrino 
mass squared differences, as described in [551 ] . 



\R\ 



' ^ m atm 
A ™sol 



The measured values adopted in [57] are given by 



(26) 



Am sol 



(7.50 ±0.20) x 10~ 5 eV 5 



Am 2 tm = (0.00232)+™eV 2 



sin 2 (2# 12 ) = 0.857 ±0.024 
sin 2 (2# 23 ) > 0.95 



sin 2 (26 



13 J 



0.098 ±0.013 



(27) 



58 



59]. 



in the 90% C. L. The global analysis for such quantities can be found in 

Suppose, for simplicity, we leave the charged lepton sector untouched and correct neutrino 
sector only. Moreover, we keep the mixings of v-i with v\$ forbidden, so that V£ is modified 
to 



I cos(f + 5) -sin(f + 5) \ 



Vr l 



(28) 



^sin(f + 6) cos(f + 6) J 
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For small 5, cos(f + 5) ~ (l/>/2)(l - 5) and sin(f + 5) ~ (1/V2)(1 + 5). From 
Vpmns = (^)V/ 





f 1 1 


1^ 




/ 


(1-5) 


-(i + <^ 


1 






1 








7! 


1 w 


u 2 






1 







v l w 2 






V 


(1 + 5) 


(1-5) , 



(29) 



we see (13) element of the PMNS matrix is given by 

25 



W, 



<3 



\/6 



(30) 



If such corrections are entirely present in the right-handed neutrino Majorana mass term 
while Y v is untouched, there would be no observable charged lepton flavor violating process. 
For example, let us introduce a doublet flavon \2- Then, its VEV modifies the diagonal 
elements of the Majorana mass matrix. With (% 2 ) = x 2 (l, 1), diagonal term has a correction 
x 2 [2N 1 Ni — N 2 N2 — N 3 N 3 }. In principle, by introducing several doublets with different VEVs, 
each diagonal term can be different. 

A. Model I 

Besides putting correction to M^r, one can find doublet VEVs giving corrections to Y v 
to make a sizable 6 Vi while dangerous charged lepton flavor violation is suppressed. To see 
this, consider the general S4 doublet VEV, (a, b) where a and b are complex numbers. With 
this VEV and coupling Ai, Y v can be modified as 

/ 



l + Ai(a + 6) 

1 + X^buj + aoo 2 ) 

1 + X^ku 2 



aw 



\ 



/ 



(31) 



In this case, YJY U in the charged lepton mass basis is given by 

1 1 + A 2 (|a| 2 + \b\ 2 ) Xx(a* + b) + X 2 ab* A x (a + b*) + X\a*b ^ 

Ai(a + b*) + X 2 a*b 1 + X 2 (\a\ 2 + \b\ 2 ) A x (a* + b) + X 2 ab* ■ (32) 
, Xt(a* + b) + X\ab* Ai(a + b*) + X 2 a*b 1 + A 2 (|a| 2 + \b\ 2 ) , 

If Ai(a* + b) + X\ab* = 0, all the off diagonal elements vanish. For example, Ai = 1 and 
a = b = u is the case. This condition also implies that off diagonal terms of Y^Y y Y^Y u vanish 
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so we do not expect any sizeable cLFV. However, this condition requires a cancellation of 
two different flavon contributions and is considered as a serious fine tuning different from 
vacuum alignment. We do not pursue this possibility any longer in this paper. 

If a* = —b and both \a\ and \b\ are smaller than one, the (12) element of Y v is given by 
— A 2 (a*) 2 . The (23) element is the same and the (13) element is its complex conjugate, — A 2 a 2 . 
In this way, LFV is suppressed quadratically even though it does not vanish. For YjY v YjY u 
term, the (12) element is 2[1 + Af(|a| 2 + |6| 2 )][Ai(a* + b) + \\ab*} + [Ai(a + &*) + X 2 a*b] 2 . 
For YJYv term, the (23) element is the same and the (13) element is its complex conjugate. 
When a* = —b, it is — 2A 2 a 2 (l + 2Af |a| 2 ) + A^(a*) 4 , which is quadratically suppressed for 
small a. For illustration, suppose \\a = X±b = ip. The stabilization of such doublet VEV is 
discussed in Appendix B. The neutrino Dirac Yukawa has the form of 



( 



V 



1 + 2ip 
1 - ip 
1 - ip 



\ 



) 



(33) 



and off-diagonal terms of YjY v in the charged lepton mass basis is suppressed to 0(p 2 ), as 
expected, 

/ l + 2p 2 - 2 - 2 



(vDHyJym 



V 



p- 



p p 

1 + 2p 2 p 2 
p 2 



(34) 



2p 2 J 



With this Y v , neutrino mass matrix is given by 

/ 



My 



\Vv\ 



, v 2 sin 2 (5 1 



2w 1 



1 + Aip- 4p 2 


-x(l + ip + 2p 2 





(l-a: 2 )(l-2zp-p 2 ) 




-x(l +ip + 2p 2 ) 


l-2ip-p 2 t 
(35) 



and the deviation of mixing from 7r/4 is given by 

-Up + 3p 2 



4x(l + ip + 2p 2 



3p 
2x 



(36) 



such that 



\V f 



e3 I 



3p 
\/6x 



(37) 
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To the first order in p, mass eigenvalues are given by 



l9 w 2 sin 2 B /l + ip 

\vA — f - — -,i 



2ip, 



ip 



(38) 



2wi V 1 + x ' 1 — x 

Taking absolute values of these eigenvalues, we obtain neutrino masses — [|?/„| 2 ?; 2 sin 2 j3/(2wi)](l/ (1- 
x),l,l/(l-x)) + 0(p 2 ). 

In summary, we expect that even though the charged lepton flavor violating effects are 
generated in the Ae term at one loop and in the m\ term at two loop, they can be suppressed 
by extra small expansion parameter p proportional to #13. With the vacuum alignment of 
the doublet flavon i(v,v), it is possible to cancel the first order correction of p and the off- 
diagonal elements of the slepton mass squared would have p 2 suppression as a result. Fig. 
[8] shows how measured #13 can be explained for the choices of p and x parameters satisfying 
observed neutrino mass squared ratio, \/R. The observed #13 ~ 0.15 can be accommodated 
for p ~ 0.1. 



B. Model II 



Of course, 8i 3 can come from both Majorana mass correction and neutrino Dirac Yukawa 
correction. Only the neutrino Dirac Yukawa coupling can affect the cLFV. To see the two- 
parameter case, consider 



l + 2ip 
1 - ip 
I 1 - ip I 



and 



/ 



M, 



N 



W\ W2 
w 1 



\ 



^ W 2 (1 - () J 



The neutrino mass is given by 



/ 



M u = -\y v 



v 2 sin 8 1 
2wi 1 — x 2 — ( 



l + Aip-C 


—x(l + ip) 



-x(l+ip) 
(1 - x 2 )(l - 2ip) - C 

1 - 2ip I 



(39) 



(40) 



(41) 
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Q, o.oo 




FIG. 8. 0i3 with respect to p and x parameters. All points in the colored region satisfy neutrino 
oscillation experiments. Neutrino #13, indicated on contour label in radian, is measured as 0.144 < 
#13 < 0.160 in la level, 0.127 < 6 13 < 0.174 in 3a level. 



where x = W2/W1 again. Then, three neutrino mass eigenvalues are given by 



2 v 2 sin 2 f3 (1 + ip ( ( 1 + ip ( 



2w x \l + x 2(1 + x) 



1 - 2ip, 



1-x 2(1 -x) 



and 



6 



W + c 2 

Ax 



Hence, we see the (13) element of the PMNS matrix is given by 

25 ^/3Qp 2 + C 2 



it; 



c3 1 



V6 



2y/6 



x 



(42) 



(43) 



(44) 



Moreover, m 2 L from Yukawa mediation is controlled by the parameter p only and (V[y (YjY u )V[ 
is the same as the previous case, 



/ 



1 + 2p 2 p 2 p 2 
p 2 1 + 2p 2 p 2 

2 J2 1 1 nJ2 



\ 



P 



(45) 



f 1 + 2?) 
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In the limit of ( — > 0, both #13 and cLFV come from the neutrino Dirac Yukawa which 
corresponds to the Model I. In the opposite limit, p — > 0, #13 is entirely obtained from 
Majorana mass term and cLFV does not appear. 

In addition, we can also constrain absolute mass scale of light neutrinos. The most 
stringent constraint on neutrino absolute mass is given by CMB data of WMAP experiment, 
combined with supernovae data and data on galaxy clustering, Ej-rrij < 0.68eV, 95% C.L. 
Conservatively, we set the bound 2.6 x 10 14 GeV < M^. Throughout paper, we use = 
5 x 10 14 GeV, the heaviest right-handed neutrino mass. 



IV. CHARGED LEPTON FLAVOR VIOLATION 



Since flavor structures of supersymmetric particles can be different from those of SM 
partners, flavor number is easily violated in SUSY. In general, the structure of the slepton 



60, 



61] . In our 



mass matrix raises dangerous cLFV. Such cLFV in SUSY is studied in 
model, once the identity structure of the neutrino Yukawa coupling Y v is broken, cLFV is 
produced. As a possible modification, one may put off-diagonal terms into Y v . On the other 
hand, when the degeneracy of Y u is broken, the combination of Y u s in the charged lepton 
mass basis, (V^)* (YjY v )V[ has off-diagonal terms as shown in Sec. IIIH The slepton mass 
squared gets extra contribution from neutrino Dirac Yukawa interactions, 



5m 2 L 



(4vr) z 



Tr^Yj] + 3Tr\Y v Y& - 3g\ - -gfjYX + 3Y^Y u YX 



(46) 



In the charged lepton mass basis, (V^)^m|V^ has off-diagonal elements and cLFV appears. 
Even though this is a general feature, it is also possible to find some parameter space in 
which charged lepton number is conserved. For example, in Sec. IIII Al off diagonal terms 
of the slepton soft mass squared, (m|)i2 can vanish for specific value of y v . Corresponding 
condition would be 



(8ml) 12 cc Tr[YX] + 3Tr[y^] - 3gj - -g\ {Y}Y v )n + 3{Y^Y V Y^Y V ) 12 = 0, (47) 



which is equivalent to 



(5m 2 L ) 12 oc 3(1 + 2p 2 )yl + 3y 2 - 3g 2 2 - -gf 

1 



- vl 



y 2 v p 2 + 3yt2p 2 



3y f 2 - 3g 2 2 - -g\ 



? + 9ytp 2 + 0(p 4 ) = 0. 



(48) 
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Near the GUT scale, g\ ~ g\ ~ 47r/28, and y t ~ 0.5 so off diagonal term vanishes for 
y v ~ 0.28. For this value of Y u , there would be no unwanted cLFV. This is different from 
the condition that diagonal contribution involving Y u vanishes, 



(6m 2 L ) u oc [3(1 + 2p 2 )y 2 u + 3y 2 - 3g 2 2 - -g 2 \yl(l + 2p 2 ) 



3^(1 + 8p 2 



3y 2 - 3gj - -g\ 



Gyt + 0(P 2 ) = 



(49) 



which is satisfied for y u ~ 0.34. 

Of course, it does not mean that y v should take the lepton number conserving value. We 
have many constraints on y v from various observations. In this paper, we try to explain the 
125GeV Higgs mass with large A term generated from y v . On the other hand, one may try 
to explain deviation of muon g — 2 from the SM prediction. Moreover, degeneracy breaking 
parameter p is used to explain sizable #13. However, it is also difficult to find an appropriate 
value of y u which satisfies all of them. In this section, we present the cLFVs for parameters 
explaining the 125GeV Higgs mass with large A term and 6 r3l . Thereafter, we visit the muon 
g — 2 constraints and the relation among # 13 , cLFV, and the Higgs mass. 



A. Experimental status 



The current experimental bounds and future sensitivities for various cLFV processes in 



the 90% C. L. are summarised in Table |H157 



62). 



B. h -> ^7 



The amplitude for L — > is written as 



T = ee^uiip - q) q 2 ^(A^P L + A?P R ) + m Xj ia , v q L \A L 2 P L + A^P R ) Uj (p). (50) 



On the mass shell (q 2 — > 0), gauge invariance imposes that the chirality preserving part does 
not contribute to the lj — > ^7 process. Hence, chirality flipping should take place in the 
on-shell lj — > la process. The decay rate is given by 



r(/^/,7) = I ^mf j (|^| 2 + |^| 2 ) 
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(51) 



Observables 
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O(10~ 13 ) [63] 


1 )1 ! 7 - 




4.4 x 1(T 8 [64] 


2.4 x 1(T 9 [69] 


R r / T _ 

1 J I I / 




3.3 x 10 _8 [64] 


3.0 x 10~ 9 [69] 


Br^ - 


■» 3e) 


l.U X 1U [OOj 




Br(r - 


■> 3e) 


9 7 v 1 fl _8 [fifil 


9 ^ v 1 n -10 IfiQl 


Br(r - 


>3/x) 
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8.2 x 10~ 10 [69] 


r(^tTi->eTi) 


4.3 x 10~ 12 [67] 


O(10- 18 )[71] 


r(^tTi— ^capture) 
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>eAu) 


7.0 x 1(T 13 [68] 




r(/jAu—> capture) 





TABLE II. Various LFV experimental bounds and future sensitivities. The table is adopted from 



and the branching ratio yields approximately 

Br(/ i -»■ la) 



a 



(52) 



In the mass insertion scheme, the chirality flipping can be easily analyzed. Consider first 
the case of the neutralino-charged slepton internal loop, as shown in Fig. [9j Fig. [9] (a) 
shows the chirality nipping from a fermion mass insertion in the external lepton line. In 
Fig. (b), chirality flipping takes place in the slepton internal line through the LR mixing 
insertion, rrij(A — /itan/3). This term consists of flavor universal part — rrij /a taxi /3, which 
can be enhanced in the limit of large tan /3 and large /i. The chirality flipping in Fig. [9] (c) 
is given by the Yukawa coupling of the lepton-slepton-Higgsino vertex. This vertex contains 
1 / cos (3 factor which combines with a sin j3 in the Higgsino-gaugino mixing to give a tan /3 
dependence to the diagram. Therefore, this diagram is enhanced in the large tan/3 limit. 
Note that it is inversely proportional to the /i, the Higgsino mass. Since this diagram contains 
SUSY mass scale only, unlike other diagrams proportional to the Higgs VEV v through rrij, 
it is dominant over all other diagrams with the neutralino-charged slepton internal loop in 
many cases. However, since the Higgsino-bino mass insertion Mz sin /3 sin 6\y and Higgsino- 
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(c) 

FIG. 9. Feynman diagrams for lj — > 1^ process with neutralino-charged slepton internal lines in 
the mass insertion scheme. 

wino mass insertion — Mz sin /3 cos 9y/ have the opposite signs, slight destructive interference 
occurs. 

Next, the case of the chargino-sneutrino internal loop is shown in Fig. [TUJ Diagrams are 
similar to those of the neutralino-charged slepton internal loop, except the absence of the 
slepton LR mixing, since the right handed neutrinos are already integrated out. Chirality 
flipping can occur either in the external lepton line (Fig. Hni(a)) or in the lepton-sneutrino- 
Higgsino vertex(Fig. [10] (b)). The latter diagram dominates over the former one, and since 
it does not have a destructive interference, it becomes the leading contribution over all other 
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a 



(b) 



FIG. 10. Feynman diagrams for lj — > Iff process with chargino-sneutrino internal lines in the mass 
insertion scheme. 



diagrams in many cases. The similar argument also applies to the discussion of muon g — 2, 
whose SUSY contribution comes from the same diagram with flavor conservation. Following 
this diagram, SUSY enhances the muon g — 2 for positive /^[nl Q]. 



In Fig. HU we show branching ratios of various lj —> Iff processes for Sec. IIII Al In the 
graph, neutrino Dirac Yukawa couplings are fixed to be y v = 0.65 and p = 0.1, while tan/3 
and SUSY breaking scale are varied. Since off-diagonal terms of m 2 L in the charged lepton 
mass basis are the same, normalized branching ratio, T(lj — > lff)/m^ are almost identical. 
Therefore, branching ratios are closely related to the total decay rate of mother particle. For 
example, since total decay rate of tau is about 5.3 times larger than that of muon, branching 
ratio of Br(/x — > ej) is about 5.3 times larger than Br(r — > e^y) and Br(r — > yrj) which are 
almost the same. 



C. r -> IT ir if 



In many cases, dominant contribution comes from the photon penguin. Z boson pen- 
guin is suppressed in general because of the accidental cancellation when the neutralino 
or chargino is pure gaugino or pure Higgsino 75[ . Such accidental cancellation is broken 
by introducing TeV scale physics which couples to the sneutrino with the sizable coupling. 
This can be realised in the R-parity violating model or in the TeV inverse seesaw, for 
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FIG. 11. Branching ratios of \i — > ery, r — > fj,j and r — > with respect to the lightest selectron 
mass for tan f3 = 3, 10, 30, y v = 0.65 and p = 0.1. 



example [72J, |75 |. 
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In our case, photon penguin is a leading contribution, so we have a simple relation between 
Br(Zj -»■ Z<7), 



Br(/j -)> ^7) 37TV mf, 4/ 

The box diagram is sup pressed in general, except for some special cases, such as in SUSY 
with Dirac gauginos 76|. 

In Fig. [T21 we show branching ratios of various lj — > 3/j processes for for Sec. IIII Al 
Fixed parameters are the same as lj — > process. We see that Br(/x — > 3e) is about 
0.018 times suppressed than Br(/i — > cy) so Eq. fl53l) is satisfied. Photon penguin is a 
leading contribution for \i — >• 3e process. In the absence of special characteristic which 
can overcome the natural size of the branching ratio, Br(7~ — >■ l^l^lf) * s c^/ 71 " suppressed 
compared to Br(/J — > / 4 r 7). 



D. [i — e conversion 

Conversion of the stopped muons in a nuclei to the electron is a promising channel to 
look for the charged lepton flavor violation. In principle there are many different operators 
including scalar, photon mediated vector, Z-boson mediated vector operators in addition to 
the dipole operator. Muon to electron conversion rate is conventionally normalised by muon 
capture rate. 

b ^ - (54) 

Here Z is the atomic number of the atom. Different target provide a different B^ e (Z) and 
relative ratio of B fl _> e (Z) of at least two different target can provide information on possible 
types of the operators as different operators predict different ratios. In supersymmetric 
modelsj?^], dominant contribution to [i — e conversion comes from the dipole operator. As 
a result, B(fi — > e)(Z) is predicted to be suppressed by a/n compared to B(fi — > ej). For 
different choice of Z, the conversion is suppressed by 10~ 3 ~5x 10~ 3 . Current limit on the 
conversion rate is comparable to // — > cy branching ratio, but the future experiments on /i 
to e conversion will have better sensitivity. We plot fi — e conversion rate with the expected 
future sensitivity of planned experiments in FigJTBl 
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FIG. 12. Branching Ratios of fi — >■ 3e, r — > 3e and r — >■ 3// with respect to the lightest selectron 
mass for tan ,9 = 3, 10, 30, y v = 0.65 and p = 0.1. 
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FIG. 13. fi — e conversion rate with respect to the lightest selectron mass for tan/3 = 3, 10,30, 
y v = 0.65 and p = 0.1. 

E. Correlation between Muon g — 2, #13, cLFV and the Higgs 

The anomalous magnetic moment of muon (muon g — 2) has a long standing sizable 
deviation from the SM prediction. The observed value is 78 1 



a M (Exp) = 11659208.9(6.3) x 10 



-10 



whereas the SM prediction 79[ is given by 



(55) 



a M (SM) = 11659182.8(4.9) x 10" 10 
so we may have new physics contribution explaining the 3.3 a discrepancy, 



(56) 



5a^ = a M (Exp) - a M (SM) = (26.1 ± 8.0) x 10 
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In the context of SUSY 



80 



(57) 



81] . muon g — 2 has the same Feynman diagram structure as 



the cLFV process /x — > ej. The crucial difference is that muon — 2 is flavor-conserving 
process, while /1 — > cy violates lepton flavor, L M and L e . Therefore, Sa^ and Br(/i — > cy) 
have a strong correlation^ 



K USY WKM 2 



Br(/i^e 7 )^3xl0- 5 (^ z ^) ( 



SUSY 



(58) 
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FIG. 14. Branching ratio of \x — > ej as a function of #13 for tan/3 = 10, y u = 0.62, p = 0.1. Future 
MEG expected bound is O(10 -13 ), we set the value 2 x 10~ 13 . Observed muon g — 2 discrepancy is 
about (2.25 ± 1) x 10 -9 , we draw Br(/i — > ej) at each muon g — 2 contribution. Green and yellow 
band indicate lcr, 3a level of neutrino #13, respectively. In upper figure, #13 is purely obtained 
from neutrino Dirac Yukawa splitting. In lower figure, only 1/15 portion of #13 is obtained from 
neutrino Dirac Yukawa. 

Moreover, the neutrino Dirac Yukawa Y u contains information on the neutrino oscillation 
observables. Since we consider the model where parameters of Y v are related to 813 and 
cLFV, we have a strong correlation between Br(/i — >• ej), # 13 , and muon g — 2 as discussed 

32 




FIG. 15. Contour plot of Higgs mass(red solid line), cLFV(black dashed line), and muon g — 2(blue 
dashed line) in Bn - y u plane for p = 0.1, tan f3 = 30. 




Fig. [TH summarises the result. Both muon anomalous magnetic moment and cLFV 
is a function of tan/3/M 2 in which M is the typical supersymmetry breaking scale. The 
cLFV has extra suppression proportional to (m|)i2. The S4 flavor model discussed here is 
constructed from the neutrino Dirac Yukawa matrix which is proportional to the identity 
matrix and does not provide any off-diagonal entry in the slepton mass squared matrix if 
#13 vanishes. Recently measured #13 ~ 0.15 provides an extra information depending on the 
origin of modification for nonzero #13. 

If the full 6 , 13 is explained by the degeneracy lift of the neutrino Dirac Yukawa matrix 
and if the entire discrepancy of the muon anomalous magnetic moment should be explained 
by light slepton, the current MEG bound tells that #13 should be smaller than 0.01 which is 
incompatible with the observation recently made. The parameter space which is consistent 
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with [i — >■ e7 bound and the #13 predicts that muon anomalous magnetic moment is at least 
1/20 times smaller than what is needed. 

If the nonzero #13 is entirely generated by modifying the neutrino Majorana mass matrix, 
there would be no cLFV even for the sizeable #13. In reality, the subleading corrections 
in the simple flavor model would appear in both sectors and the observed #13 would be a 
combined result of various sources. The bottom plot of Fig. [14] shows the hybrid case in 
which only 1/15 of the #13 is from the neutrino Dirac Yukawa modification. In this case #13, 
muon anomalous magnetic moment can be explained at the same time. The \i — > bound is 
satisfied and the the consistent region can be reached by the planned future MEG experiment 
as it predicts larger branching ratio of \i — > ej than the planned expected sensitivity. 

Fig. [J5] shows the tension between the muon g — 2 and the Higgs mass. Even if we take 
the model in which the neutrino Dirac Yukawa remains to be proportional to the identiy 
matrix such that no cLFV constraints apply, 125 GeV Higgs mass needs Bn much larger 
thna 300 TeV. Then the slepton is too heavy and the muon g — 2 is much smaller than 10~ 9 . 
The figure also shows an interesting feature that the off-diagonal elements of m\ vanish at 
around y v = 0.3 and the cLFV bounds are very weak at around y v = 0.3. 

V. CONCLUSION 

We considered the right-handed neutrinos as messengers of supersymmetry breaking in 
the minimal gauge mediation. Direct coupling of neutrino messenger with the Higgs field 
H u and the lepton doublets L, provides soft-trilinear A term for the top Yukawa and can 
help increase the light Higgs mass by realising the maximal stop mixing scenario. We call 
this setup as 'neutrino assisted gauge mediation'. The Yukawa mediation given by neutrino 
messengers also appear at soft scalar masses of the Higgs H u , the lepton doublets L,;. At 
the same time it affects the soft scalar masses of the fields which couple to H u and L, at two 
loop. Among those, the stop mass squared gets the largest correction as the top Yukawa 
coupling is of order one. For y v slightly larger than 0.7, the correction is big enough to make 
stop tachyonic. Therefore, this realises the natural supersymmetry spectrum. At the same 
time maximal mixing is achieved by two effects, large A t and small m~ at around y u ~ 0.7. 
In general this effect allows to explain the observed Higgs mass at around 125 GeV using 
around 1 TeV stop mass. Compared to the case when the neutrino assistance is turned off 
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(y v = 0), about 5 GeV of the Higgs mass is enhanced. 

In general the off-diagonal entry of the slepton mass squared, m|, appears at the mes- 
senger scale and can make the charged lepton flavor violating process to occur. The detailed 
quantitative prediction of cLFV highly depends on flavor model building. We provided a 
representative model based on S4 flavor symmetry in which the Dirac neutrino Yukawa can 
be proportional to the identity if 6* 13 = 0. For nonzero #13, two options were considered. 
Firstly, the total 6 13 can be explained by the modification of the neutrino Dirac Yukawa 
matrix. Secondly, the #13 can be explained by modifying the Majorana mass matrix of neu- 
trinos. For the former, very stringent bound on the slepton mass comes from /i — > e^y and 
the slepton should be heavier than 2 ~ 4 TeV, depending on tan 0. Also for the slepton 
mass at around 2 TeV with tan /3 = 10, the /1 — > is just below the current experimental 
bound and we expect to observe the fi — > ej in the near future. 

Even for the second case in which we can safely avoid cLFV constraints, the neutrino 
assisted gauge mediation (in its minimal form with one copy of 5 and 5 messenger) sets the 
lower bound on the slepton mass to explain the Higgs mass. 1 ~ 2 TeV slepton mass at the 
same time sets an upper bound on the possible contribution to muon anomalous magnetic 
moment and a M ~ 10" 10 is the upper bound. 

In this paper we proposed the neutrino assisted gauge mediation and showed a possible 
way out to avoid the strong cLFV constraints. Even then the current scheme has a tension 
with the muon anomalous magnetic moment which needs a lighter slepton. The extension of 
the minimal neutrino assisted gauge mediation to multiple messengers might ameliorate the 
tension between the spectrum needed to explain the Higgs mass and the muon anomalous 
magnetic moment. 
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(tan p = 10, = 360 TeV) 


(tanp = 30, B N = 300 TeV) 




2957, 2961, 3013 


2429, 2465, 2502 


ei,Ati,fi 


1364, 1364, 1333 


H39, H38, 880 




3013, 2962, 2954 


2503, 2467, 2427 


Ul,Cl,tl 


2827, 2827, 634 


2384, 2384, 637 


di,h,h 


2853, 2853, 2820 


2406, 2406, 2283 


U2,C 2 ,t2 


3177, 3177, 2252 


2675, 2675, 1868 


d2,S2,fo 


3178, 3178, 2297 


2676, 2676, 1893 


ho,A,H ,H± 


125, 1705, 1705, 1707 


125, 1031, 1031, 1034 


Xl,X2,X3,X4 


487, 850, -892, 980 


405, 713, -758, 829 


X+,X- 


849, 980 


712, 829 


a 


2514 


2126 



TABLE III. Sparticle spectrum at the point giving 125 GeV Higgs mass with the lowest Bn 
APPENDIX 0: SPARTICLE SPECTRUM SAMPLE POINT 

APPENDIX A: REPRESENTATIONS OF 5 4 SYMMETRY AND TENSOR PROD- 
UCTS 

5*4 is a non-abelian discrete symmetry and consists of all permutations among four quan- 
tities. For a review, see jssj]. Irreducible representations of 5*4 are two singlets 1,1', one 
singlet 2, and two triplets 3, 3'. Tensor products among them are given as follows: 



^ xi \ ( yi ' 

X 2 

\ X3 /s 



in 



(Xiy! + X 2 V2 + ^32/3)l + 



Xxyx + UJX2V2 + (*rx 3 y 3 
xiyi + uj 2 x 2 y 2 + uxm 



( x 2 V3 + x 3 y 2 




^ x 2 yz 


- x 3 y 2 ^ 


xsvi + xm 


+ 


xm 


- xm 


l xm + x 2 yi j 


3 




-x 2 yi j 



(59) 
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(tan /? = 10, = 240 TeV) 


(tan p = 30, B N = 200 TeV) 




1971, 1974, 2009 


1657, 1682, 1707 


ei,Ati,fi 


915, 915, 894 


770, 769, 590 




2010, 1976, 1970 


1709, 1684, 1657 


Ul,Cl,tl 


1937, 1937, 521 


1633, 1633, 404 


di,h,h 


1954, 1954, 1931 


1650, 1650, 1564 


U2,C 2 ,t2 


2169, 2169, 1569 


1828, 1828, 1286 


d2,S2,fo 


2170, 2170, 1586 


1829, 1829, 1291 


ho,A,H ,H± 


123, 1220, 1220, 1223 


123, 679, 679, 684 


Xl,X2,X3,X4 


322, 600, -729, 757 


267, 466, -520, 579 


X+,X- 


600, 757 


465, 579 


a 


1737 


1470 



TABLE IV. Sparticle spectrum at the point giving 123 GeV Higgs mass with the lowest Bn 







H 


2 2 


X 


2/2 


[x 3 j 


3' 


{y 3 ) 



(SiJ/i +222/2 +£32/3)l 



ZlJ/l + WX2J/2 + W 2 X 3 J/3 

xij/i + w 2 x 2 j/2 + ujx 3 y 3 



222/3 + 2 3 2/2 
2 3 2/l + 212/3 
^2i£/2 + 222/1 y 



^ 2 2 l/ 3 - 2 3 2/2 ^ 



(60) 



232/1 - 2lJ/3 

2i?/2 - x 2 yi 



J 3' 



' ^i I I 2/1 ' 

22 

\ X3 J 



x 



2/2 



(xiyi + 222/2 + 232/3)1/ + 



xiyi + CJ222/2 + w x 3 y 3 
-(xxyx + w 2 x 2 2/ 2 + WI32/3) 



+ 



2 2 2/3 + 2 3 2/2 
2 3 2/l + 2l2/3 
2l2/2 + 2 2 2/l 



' 2 2 2/3 - 2 3 2/2 ^ 



(61) 
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2 3 2/l - 2l2/3 
2l2/2 - 2 2 2/l 











H 











(xiy 2 + x 2 yi)i + {xiy 2 - x 2 yi)v + 




(62) 




2/2 



(xi + x 2 )yi 
(u 2 xi + ux 2 )y 2 
y {uxi + u 2 x 2 )y 3 J 



(xi - x 2 )yi 
(u 2 xi - ux 2 )y 2 
y {uxi - u 2 x 2 )y 3 J 



(63) 




x 



1)2 
V 3 J 



(xi + x 2 )yi 
(cu 2 xi + ux 2 )y 2 
y (ujx 1 + u 2 x 2 )y 3 J f 



(xi - x 2 )yi 
{u 2 xi - ux 2 )y 2 
y {uxi - u 2 x 2 )y 3 J 



and trivially, we have 3xl' = 3', 3' x V = 3', and 2 x V = 2. 



(64) 



APPENDIX B: REMARKS ON THE FLAVON VACUUM STABILITY 



In 55], it was shown that v4 4 triplet flavon vacuum in the direction of (1,1,1) and 
(1,0,0), ((0,1,0), (0,0,1) are the same) is favored compared to other directions, such as 
(1,1,0). Since symmetry is the subgroup of the S4 composed of even permutations, 
similar arguments hold. In this Appendix, we argue that triplet flavon directions favored 
in A4 model are also favored in the S4 model and that S4 doublet vacuum favors the (1,1) 
direction. 

Rigid SUSY makes the discussion more simple, because the potential V is minimized at 
(V) = 0. On the other hand, extra symmetries like Z4 and U(1)l more restrict possible 
terms in the superpotential. Suppose that U(1)l symmetry is discretized to, for example, 
Z 8 symmetry. In this case, only quartic terms <3> 4 and x 4 are allowed. Let us assume that 
breaking of extra symmetries introduces quadratic term, like m^ 2 or m 2 x 2 - To achieve 
this, let us consider 'Z 4 breaking singlets' ipi, z/^ and 'lepton number breaking singlets' ip 2 , 
ip 2 with 1S4 x Z4 x U(1) L quantum numbers 



^: (1,3,0), (1,1,0), 
^: (1,0,2), ^2: (1,0,6). 



(65) 
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They do not combine with ELHd, NLH U and NN to make singlets under all symmetries 
imposed. Then, they can couple to $ 2 or x 2 such that a superpotential is given by 

Wfa, ^1,^2,^2) =^-[$ 2 ^iVi + x 2 ^H 

- M 1 ^ 1 + ~[«i$iVi) 2 + ^if + ^i) 4 ] (66) 

In this superpotential, ViV'i an d ^2^2 pahs have VEVs and they provide m^ 2 + m 2 x 2 
terms. With this setup, the triplet superpotential has the form of 

W = mS 2 + K x 2 +y 2 + 2 2 )2 + A2 (x 2 + ^2 + w 2^2 )(x 2 + ^2 + ^2) 

, A A (67) 

+ ~J^\ X V + V z + ^ x ) 

where 5 = (x, 2;) represents the generic S4 triplet such as $ or x- Note also that the 
superpotential has an accidental Z 2 symmetry under which and 7/^2 are odd whereas 
other fields are even. If this Z 2 symmetry is imposed, (§ 2 ipi/A 3 )ELH<i and (& 2 ip 2 /A 2 )NN 
terms, which change the flavor structure in the subleading orders are forbidden. In this 
case, charged lepton Yukawa coupling structure in dimension-4 operator is preserved up to 
dimension-6 operator whereas Majorana mass structure in dimension-3 operator is preserved 
up to dimension-5 operator so corrections to them are highly suppressed. 

Each term of the F-term potential V = \dW/dx\ 2 + \dW/dy\ 2 + \dW/dz\ 2 is given by 

dW 4Ai , 2 2 2\ 2A2 , 2 2\ 2A3 . . . . 

— — = mx H — — x[x + y + z ) H — — x(2x —y — z ) H — — (y + z)[xy + yz + zx) 
ox A A A 

dW 4Ai 2 2 2\ 2A2 2 2 2\ 2A3 . 

— — = my H — — xur + ?/ + z ) H — —y(2y —z — 2 H — — (z + x)(xy + yz + zx) (68) 
oy A A A 

dW 4Ai , 2 2 2\ 2A2 /„ 2 2 2\ 2A3 , w . 

— — = mz + — z(i +y + z ) + —^-z{2z —x — y ) + —^{x + y)(xy + yz + zx). 

Stable vacuum requires that these three terms should be zero simultaneously. For vacuum 
(S) = v(l, 1, 1), three terms give the same condition, 

3 

12(Ai + A 3 )Q-) +mv = (69) 

so the vacuum is stabilized at v 2 = — mA/[12(Ai + A3)]. For vacuum (S) = v(l,0, 0), the 
second and third terms vanish trivially and the first term gives 

3 

4(Ai + A 3 )Q-)+7m; = (70) 
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so the vacuum is stabilized at v 2 = — mA/[4(Ai + A3)]. The vacuum in the direction (0, 1, 0) 
and (0, 0, 1) gives the same result by permutational property of 5*4. On the other hand, 
vacuum (S) = v(l, 1, 0) gives two conditions, 

^-(8Ai + 2A 2 + 2A 3 )+7m; = 

A (71) 
\ 3 v 3 = 0. 

If A3 is not forbidden by another symmetry, v — is the only solution and nontrivial vacuum 
can not be developed. 

S4 doublet stabilization can be discused in the same way. Renormalizable superpotential 
for doublet (x, y) is written as 

W = m(xy) + A(x 3 + y 3 ) (72) 



and stabilization condition 



^ = 2my + 3Xx 2 = 
2mx + 3Xy 2 = 



dy 

requires that x = y. So the vacuum choice for Eq. is stable 



APPENDIX C: COMMENT ON KAHLER POTENTIAL CORRECTIONS 

In our setup, Yukawa couplings are constructed from non-renormalizable dimension-4 
superpotential with several flavons. These flavons also appear in the non-renormalizable 
Kahler potential and kinetic terms are written in the form of 



(74) 



where <ft and ip represent bosonic and fermionic fields, respectively. The Kahler potential of 
charged lepton supermultiplet L is given by 



K 



$t$ 

1 + at— -r- + a 2 
A z 



A 2 



S4 singlets 



+ ■ 



(75) 



and similar terms can be written for other fields, E^E, N^N, W u d H u ^, etc. Then we have 



quite complicate terms. For example, from ($g$ 3 /A 2 )L^L where $3 vacuum is given by 
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vz(l, 1,1), we have 

54 singlets 



A 2 



,i~^(L\Li + L2L2 + L\Lz) 

(76) 



Since ($)/A = U2/A is responsible for charged lepton Yukawa couplings, we see 47rf 2 /A > 
Y T = m T / \{v / \/2) cos/3] ~ 0.1 for tan/3 = 10. On the other hand, Xs nas another vacuum 
direction, w 2 (0, 1,0). Then 



,2 



— a 2,i -r| {L[L\ + + -^3-^3) 

54 singlets A (77) 

+ a 2,2T-|( — L\Li + -^2-^2 — L\L s 



so it just rescales the fields. Moreover, since See-Saw scale is about 10 14 GeV, we have 
suppressed effect, 47r;\;/A ~ 0.01 with A is the GUT scale. In the same way, doublet and 
singlet flavons in the Kahler potential just contribute to the field rescalings. 

Physical fields are defined with canonical kinetic terms, so we should make field redef- 
initions and they affect flavor structures in principle. In our work, however, such effects 
are not considered by assuming small coefficients 0^2- For example, diagonalization of Ye 
demonstrated above is not affected if ai^/A 2 ) < (m e /m r ) ~ 3 x 10~ 4 , i.e. a\ < 3. 

On the other hand, mixings in the Kahler potential between flavons can be dangerous. 
For example, kinetic mixing between flavons such as iplipl^sXs/-^ 2 can introduce small 
correction to Y E or M N with unwanted S4 triplet vacuum direction. Such effect is suppressed 
by iplipl/A 2 and can be more suppressed with tiny coefficient. 
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